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Many promising schemes for quantum information processing (QIP) rely on few-photon inter- 
ference effects. In these proposals, the photons are treated as being indistinguishable particles. 
However, single photon sources are typically subject to variation from device to device. Thus the 
photons emitted from different sources will not be perfectly identical, and there will be some varia- 
tion in their frequencies. Here, we analyse the effect of this frequency mismatch on QIP schemes. As 
examples, we consider the distributed QIP protocol proposed by Barrett and Kok [3| , and Hong-Ou- 
Mandel interference which lies at the heart of many linear optical schemes for quantum computing 
010)01 Si- In the distributed QIP protocol, we find that the fidelity of entangled qubit states de- 



QQ ' pends crucially on the time resolution of single photon detectors. In particular. There is no reduction 

, in the fidelity when an ideal detector model is assumed, while reduced fidelities may be encountered 

■ when using realistic detectors with a finite response time. We obtain similar results in the case of 

' Hong-Ou-Mandel interference - with perfect detectors, a modified version of quantum interference 

, is seen, and the visibility of the interference pattern is reduced as the detector time resolution is 

^ ■ reduced. Our findings indicate that problems due to frequency mismatch can be overcome, provided 

' sufficiently fast detectors are available. 

^ : 



qubits. As a result there is now growing interest from ex- 
perimental groups in implementing distributed schemes 
[Tsl [l^ . Other applications of single photon interference 
in QIP have also been proposed, such as quantum re- 
peaters [itI [isj , and are currently being actively pursued 
by experimental groups [19j. 

Few-photon interference effects are often said to re- 
quire identical photons, such that the frequency, polar- 
ization, and temporal envelope of each photon should 
be indistinguishable. This is because these experiments 
make use of a beam splitter to erase 'which path' in- 
formation, so that the source of each photon cannot be 
inferred from the detector signal, even in principle. Any 
additional information carried by the photon (such as the 
frequency) could, in principle, be used to infer the path 
that the photon took, and therefore will tend to degrade 
the interference. From this perspective, one expects that 
few-photon interference cannot be observed between pho- 
tons from sources of different frequency. Indeed, if one 
restricts ones attention to the total coincidence rates, 
this is indeed what is observed in a Hong-Ou-Mandel 
type experiments. For sufficiently detuned single pho- 
ton sources, the photons behave as independent particles, 
each exiting either port of the beam splitter with prob- 
ability i, and no interference is observed. However, this 
begs the question, "where does the interference 'go'?" 
Usually we only expect quantum effects to vanish in the 
presence of some kind of noise or decoherence process. 

Some insight into this issue has been provided in a se- 
ries of intriguing theoretical and experimental results by 
Legero and co-workers, concerning the Hong-Ou-Mandel 
effect with different frequency photons [j, They 
showed analytically that if one can perform time resolved 
measurements in a Hong-Ou-Mandel type experiment 
with detuned photons, a type of quantum interference 
can still be observed. The probability of both photons 
appearing at the same output port of the beam splitter is 



I. INTRODUCTION 

Few-photon interference effects are of fundamental in- 
■ terest, as they have no classical analogue, and demon- 
strate the quantum nature of the radiation field. The 
archetypal example of such an effect is two-photon inter- 
ference as observed by Hong, Ou and Mandel in 1987 Q. 
There two identical photons, each incident on a separate 
input port of a beam splitter, coalesce such that both 
photons are always detected at the same output mode of 
QQ ' the beam splitter. More recently, single-photon interfer- 
psj . ence effects have been proposed as a means to entangle 
' remote pairs of atomic systems 0, Here, the atoms 
, emit photons which are then incident on a beam split- 
' ter, followed by measurements on the output ports of the 
[ — \ beam splitter. The role of the beam splitter is to coher- 
' ently erase 'which path' information. Since the photons 
^ ' are identical and the observer cannot know which atom 
the photon was emitted from, the result is to prepare the 
pair of atoms in an entangled state. 

d ' Such interference effects are currently the subject of 
much interest, as they potentially provide a route to 
scalable quantum information processing. Few-photon 
interference lies at the heart of schemes for linear op- 
tics quantum computing 0, [^, [13, [HI- This is also 
central to many hybrid light-matter quantum comput- 
ing schemes [ll, [13] , in which remote matter qubit 
systems (such as trapped atoms, quantum dots, or impu- 
rity centers in solids) can be entangled via single photon 
interference effects, in a way such that efficient quan- 
tum computation is possible. This approach can signifi- 
cantly simplify scaling the computer to a large number of 
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now no longer a constant value, but rather oscillates as a 
function of the time between the photon detection events. 
The frequency of this oscillation is given by the detun- 
ing between the two photons, and thus has been called 
a 'quantum beat' of two photons ^Sj. This effect was 
observed experimentally using successive photons from a 
single photon source implemented by a Raman transition 
in an atom-cavity system Q. 

Aside from the conceptual interest in these effects, the 
issue of interference between photons from non-identical 
sources is now of significant practical importance, since 
there is much interest in using such effects in quantum 
information processing. Scaling these proposals will re- 
quire observing quantum interference between photons 
from many different sources. These sources may be man- 
ufactured devices, such as quantum dots or other sys- 
tems coupled to micro-cavities [l^l, and as such will be 
subject to fabrication imperfections. In particular, some 
variation of the relevant optical frequencies of the de- 
vices is to be expected. This will also be a problem in 
certain 'natural' systems, such as nitrogen-vacancy cen- 
ters in diamond, which typically experience a spread in 
their resonance frequencies due to local strain fields [1^ . 
While it could be possible to tune such systems over a 
certain frequency range, it may nevertheless be difhcult 
to bring all of the sources into resonance with each other. 
Thus it is important to understand the extent to which 
mutual detuning of the sources induces errors in QIP 
schemes making use of single photon interference effects, 
and what can be done to mitigate such errors. 

In this paper we investigate the practical and funda- 
mental aspects of the effect of detuning on few-photon 
interference by considering two particular cases. Firstly, 
we consider the entangling operation introduced by Bar- 
rett and Kok This operation allows the preparation 
of entangled states of remote qubits, and furthermore can 
be used to generate graph states of multiple qubits, and 
hence is a resource for scalable, universal quantum com- 
putation. In addition we consider the Hong-Ou-Mandel 
effect. This is of interest because it is one of the best 
known few photon interference effects, and is also central 
to multiple schemes for linear optics quantum computa- 
tion 0, 13 and quantum repeaters (isj . In both cases, we 
first consider the corresponding effects with detuned pho- 
tons in the case of idealized (i.e. perfect time resolution) 
detectors, and find that a modified version of the entan- 
glement/interference effect persists. We then consider 
the opposite case, where the detectors have very bad time 
resolution, and find that in this limit, the interference (or 
entanglement) is indeed reduced. Loosely speaking, the 
degree of reduction of entanglement /interference depends 
on the extent to which the photons are distinguishable 
in the frequency domain. 

By making use of an explicit model of the photode- 
tectors [1^ [2^, we also consider the intermediate case, 
where some time resolution present in the detector out- 
puts, and quantify how the entanglement/interference is 
modified as the detector resolution varies. Our results 



are of direct practical benefit for the implementation of 
these schemes, as they allow one to determine the level of 
error that can be expected for a given detuning and de- 
tector resolution. Furthermore, we hope the results will 
aid in understanding the nature of single photon interfer- 
ence effects with detuned photons, and give some insight 
into where the entanglement/interference 'goes'. 

Although we focus on two particular examples of few- 
photon interference in this paper, the techniques are rea- 
sonably generic and can therefore also be applied to many 
other schemes which involve similar effects. The effect of 
frequency mismatched photons has also been examined 
in slightly different contexts, such as sources of entangled 
photon pairs [13] ■ We also note that a potential solution 
to the frequency mismatch problem has been proposed 
[2^ . This scheme makes use of acousto-optic modulators 
as 'frequency beam splitters' which can be used to erase 
the frequency information of the photons. 

The remainder of this paper is organized as follows. 
In SectionJIIl we review the entangling scheme of Barrett 
and Kok [l| . We examine the effect of photon detuning 
in this scheme in Section IIIII The case of ideal time 
resolution detectors is considered in Section fill Ai while 
in IIIIBI we examine the opposite limit of very bad time 
resolution detectors, and in Section llll CI the intermediate 
case is analysed. In Section llVI we explore the influence 
of detuning on the Hong-Ou-Mandel effect, again in the 
regimes of good, bad, and intermediate time resolution 
detectors. Finally, we summarize our findings and draw 
some conclusions in Section IVl 



II. ENTANGLING ATOMS 

In this section, we review the method proposed in Ref. 
[l[ for entangling remote pairs of qubit systems, which 
could be trapped atoms, ions, quantum dots, or impu- 
rity centers in solids. This method actually implements 
a non-deterministic parity measurement such that, when 
successful, a projection of the form |01)(01| -I- |10)(01| is 
performed on the joint state of the qubits. A positive 
outcome is heralded by a particular sequence of detec- 
tor clicks. If these are not observed, the operation has 
failed, and the qubits can be reset and the operation reat- 
tempted. This operation, combined with single qubit ro- 
tations and measurements, is sufficient to efficiently gen- 
erate arbitrary graph states of multiple qubits, which in 
turn permit universal measurement-based quantum com- 
putation. A number of other schemes have also been 
proposed for remote entanglement via few-photon inter- 
ference 0, i [13, [11 [3 m, [13, S m S iH, p , and 

although we do not consider these schemes explicitly, the 
results in this paper are expected to also be applicable 
to those schemes. 

The setup proposed consists of two atoms [il] inside 
separate cavities, as shown in Figure [1] which are as- 
sumed to have equal resonant frequencies lui = u>2- The 
qubit levels |0) and |1) are long lived, low- lying states of 
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FIG. 1: Schematic diagram of the setup proposed, (a) Level 
structure of the atomic system. |0) and |1) are low lying, 
long lived states representing the qubit degree of freedom. 
|1) is connected to |2) by an optical transition which may be 
addressed by a 7r-pulse of a laser to swap population between 
these states, and coupled to the cavity mode with coupling 
constant g. (b) Setup for remote entanglement. The light 
from each atom is collected (possibly with the aid of an optical 
cavity on resonance with the 1-^2 transition) and mixed on 
a 50:50 beam splitter. Photon counting detectors D+ and D- 
monitor the output modes of the beam splitter. In the original 
proposal, the frequencies u>i and u>2 of the optical transitions 
are assumed to be equal In general, this will not be the 
case and we examine this scenario in Section Hill 

the atoms. In addition there is an excited level, |2), such 
that an optical transition between |1) and |2) couples 
resonantly to the cavity mode of the respective cavity. 
The use of cavities is not strictly necessary for the ideas 
presented here, since the same protocol will work if the 
light emitted from the optical transition is monitored in 
free space. However, coupling via a cavity may offer an 
increase in the success rate of the protocol with respect 
to the same setup without cavities. The protocol for en- 
tangled pair generation proceeds as follows: 

1. Prepare atoms in the product state |+) (E) |+), 

2. apply a 7r-pulse on the 2^1 transition to prepare 
the atoms in i(|00) + |02) + |20) + |22)), 

3. monitor cavity output for a time significantly longer 
than the decay time of the 2^1 transition; a click 
in either detector signals a successful first round. 
The absence of a click implies a failure of the op- 
eration, and the protocol should start again from 
step (1). 

4. Apply a bit-flip, cTj., on the qubit states to perform 
|0) ^ |1) on each atom. 

5. Repeat steps (2) and (3); a second click in either de- 
tector signals the successful completion of the pro- 
tocol. 

We now review these steps in more detail. The atoms 
are initially individually prepared in the |-|-) = "^(|0) -I- 
|1)) state, such that the combined state of the atoms may 



be written as 

IV'(O)) = 1+) ® 1+) ^ ^(|00) + |01> + |10) + 111)) . (1) 

Then a 7r-pulse is applied to excite the 1-2 transition 
such that the state of the atoms becomes ^(lOO) -I- |02) + 
1 20) -I- 1 22)). Next the cavity leakage is monitored using 
the detectors D+ and D_. A click in either detector 
signals a successful first round of the protocol. If the 
detection process were perfect it would be possible to stop 
the protocol here and be confident of having correctly 
performed an entangling operation on the atoms. The 
final state of the atoms in this case is given by 

|^±) = -i=(|01)±|10)), (2) 

where the sign between the terms is determined by the 
detector where the click was observed. This state is ob- 
tained because the only parts of the initial state that 
create exactly one cavity excitation are initial atomic 
states |01) and |10). The presence of the beam split- 
ter erases the which-path information such that a click 
in either detector will not reveal any information about 
which atom/cavity the excitation originated from. 

However, in general the photon collection process as 
well as the detectors themselves will not be perfect. 
Therefore the detection of a single click will lead to a 
mixed state over the one- and two-excitation parts of the 
atomic states, 

p± = Pi|^±)(V'±| + (l-pi)|ll)(ll|, (3) 

where pi is the probability of there having been only 
one photon. Note that strictly speaking this neglects the 
presence of dark counts in the detectors, which would 
lead to a |00)(00| contribution to the state. However, for 
existing detectors the dark count rate is typically much 
smaller than the atom/cavity emission rates involved in 
the protocol, thereby justifying this approximation. In 
addition, detector dead-times do not affect this proto- 
col, since for a successful outcome, only a single click is 
observed on each round, and two-excitation events are 
rejected as a result of the post selection process (as we 
describe below). 

The solution proposed to overcome the presence of the 
two-excitation component in the state, was to apply a bit- 
flip pulse to each qubit, i.e. a ax operation on |0) and |1), 
followed by a second round of the protocol. The bit-flip 
operation has no effect on \tp^), but changes |11) |00). 
Therefore a second round of the protocol resulting in a 
second click eliminates the |00) part of the state as |00) 
does not couple to the relevant optical fields. Thus no 
photon can result from this component of the state. Then 
the final state of the system, conditional on observing a 
photon in each round, is given by 

IV) = 4H(ioi) + (-iriio))' (4) 
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where m = if both chcks occur in the same detector, 
or m = 1 if they are in different detectors. This state is 
maximally entangled and independent of click times. The 
independence on click times is due to the assumption that 
uji = L02- We shall consider the case of detuned cavities 
in the following section. 

The success rate of the protocol is 



1 



(5) 



where 77 is the combined efficiency of collection and detec- 
tion of the photons, while the factor ^ is from the pop- 
ulation of the initial state in the {|01), |10)} subspace. 
In spite of the inherent non-determinism of this proto- 
col, efficient quantum computing is still possible using 
this operation, in principle with any success probability 
larger than zero [1]. (The price to be paid for this is in an 
overhead cost for building the cluster states, which may 
become impractical for very small success probabilities 

Note that the fidelity of the entangling operation is not 
affected by photon collection or detection efficiency, since 
only outcomes in which a photon was observed on each 
round of the protocol are retained (i.e. post-selected) as 
successful outcomes. This is a useful fact for analysing 
the fidelity of the scheme in the presence of other imper- 
fections, since it means that one can essentially ignore 
the detection and collection inefficiencies in calculations, 
and still arrive at reliable values for the fidelity. One 
should still be mindful, however, that photon loss will 
lower the success probability. This is an approach that 
we will adopt in the remainder of this paper. 



III. ENTANGLING OPERATIONS WITH 
DETUNED SOURCES 

In this section we examine the effect of cavity frequency 
mismatch on entangling operations. In particular we ex- 
amine the scheme outlined in Section [TTl The setup is 
generalized to cavities with different frequencies, loi and 
UJ2. The atomic transition frequencies are each still as- 
sumed to be resonant with the corresponding cavity tran- 
sition, such that Wi2;i — toi and lji2-2 — ^2, where uJi2;j 
is the frequency of the 1 — > 2 transition in atom j. We 
make this assumption partly to simplify the analysis so as 
to concentrate specifically on the effect of frequency mis- 
match between two sources, but it is also reasonably well 
motivated physically. This is because in systems where 
the atomic transition frequencies are not naturally on 
resonance with the cavity resonance, it may still be pos- 
sible to tune the transitions into resonance, for example 
through the use of Stark or Zeeman shifts of the atomic 
levels. However, tuning the cavities so that they are also 
mutually on resonance may be more difficult, especially 
in the case of monolithic micro-cavities. We also restrict 
our attention to the case where both cavities have the 
same decay rate ki = K2 = k, and the same atom-cavity 



coupling strength, gi = g2 = 9- A treatment when this 
is not the case has already been presented in fl|. 

Using the Quantum Jump (Q J) formalism [3J, HI, [H, 
IstI . Issj which is particularly well suited to describing 
systems under continuous observation, we find that the 
Hamiltonian describing the evolution of the system con- 
dition on no-photon emissions is given by 



+hg{b,\2),,{l\ + 6]|1>,,(2|) - '-^b% , (6) 



where the index j — 1,2 labels the respective atom-cavity 
systems. We have defined the energy of the degener- 
ate ground states |0)j, \l)j to be zero. The first term 
corresponds to the energies of the atoms and the cavity 
fields respectively. The second term describes the Jaynes- 
Cummings interaction between the cavities and the 1-2 
transitions of the atoms, while the last term comes from 
the QJ description of the cavity-free field interaction. 
This term is non-Hermitian and leads to a decrease in the 
norm of the state vector which quantifies the decrease in 
probability of the system not emitting photons. We now 
transform to an interaction picture via the unitary evo- 
lution operator Uq — exp[-ji^i{J2'^j=i |2)jj(2| + b''^bj)t]. 
The state transforms according to \ilJi{t)) = f/p |'!/'s(i)) 
where \ips{t)) is the state in the Schrodinger picture. The 
dynamics in the interaction picture is given by i?condj 

where i?cond 
transformation to 



ihUf^Uf). Applying this 



gives 



i/,„„d = E Kgjh\2),, (1| + g*b]\l),, (2|) 



3=1 

+SA(|2)22(2|+452)-z|k6]6,, (7) 

where we have defined A = ljJ2 — uji. The associated jump 
operators which describe the evolution of the system in 
the event of an emission of a photon out of either cavity 
are given by 



Ri = Vk^i : 
R2 = \^b2 ■ 



(8) 



Now the unconditional master equation for this dissipa- 
tive system may be written as 

P ^ j- (i?condP - pHl^J + R,pR\ + R2PRI . (9) 

The effect of a 50-50 beam splitter which mixes the cavity 
outputs bi and 62, as described in Section [ll] is described 
by the beam splitter transformation, 

61 -f 62 



V2 ' 
hi - b2 
V2 ' 



(10) 
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where c+ and c_ are the two output modes of the beam 
sphtter. Note that the transformed operators, c+ and 
c_, do not have an exphcit time dependence in the in- 
teraction picture. This is a consequence of the form of 
Uq - both cavity jump operators, 6i_2, receive the same 
time-dependent phase shift, and so there is no relative 
phase between the bi and 62 terms in the transforma- 
tion. This transformation leaves the master equation un- 
changed, but will however influence single trajectories. 
This is reflected by the change in jump operators. 



R+ = + ^2) = Vkc+ 



(11) 



Now we use the fact that for the over damped (i.e. Pur- 
cell) regime, i.e. when k ^ g, we may eliminate the pop- 
ulated cavity mode to simplify the analysis [s^. This 
approximation is possible as the population in the cav- 
ity mode remains negligible in this regime. Then we find 
that the effective Hamiltonian of the system is then given 
by 

2 

J/eond = SA|2)22(2|-;^zf«;eff|2),,(2|, (12) 

J = l 

where KeS = 4g^/K is the effective decay rate of the 
atoms, for decay via the cavity mode. Similarly the ef- 
fective jump operators are given by 



R+ = \/^(|l)ll(2| + |l)22(2|), 



R- - ^^(|l)ii(2|-|l)22(2|) . (13) 

Note that Eqs. p2HT3|) are also applicable to setups with 
no cavities, where the spontaneously emitted light is di- 
rectly collected via a system of lenses and other optical 
elements [lEIiBl- In this case, KcS should be replaced by 
the appropriate rate for spontaneous emission into the 
collected mode. 



A. Ideal detector case 

In this section we consider the case of remote entan- 
gling operations in the case where the time resolution of 
the detectors, i^, is 'ideal', in the sense that it is very 
much shorter than 1/A. The time resolution can be un- 
derstood as the uncertainty in the time at which the pho- 
ton caused a change in the detector, due to technical im- 
perfections in the detector, and finite bandwidth of the 
associated electronics. Note that we cannot assume truly 
infinitesimal time resolution, since the quantum jumps 
formalism that we apply here uses Born, Markov, and 
rotating wave approximations which break down on very 



short timescales on the order of the inverse of the optical 
frequency. Thus the results in this section are valid in 
the regime l/wi^2 ^ tr ^ 1/A. 

As noted at the end of Section [TTl we can analyze the 
scheme as if the collection and detection efhciency was 
perfect, since we are only interested in the final state in 
the case when detector clicks where actually observed on 
both rounds of the protocol. Finite collection and detec- 
tion efficiencies will not affect this conditional state, but 
will simply reduce the overall success probability. In this 
case, the evolution of the system when no detector clicks 
are observed is described by the Shrodinger equation. 



-rHr, 



mi)) 



(14) 



Since ifcond is non-Hermitian, this evolution is non- 
unitary, and thus \'4^{t)) is unnormalized. The norm- 
squared of the wavefunction, ||'0(t))p, can be interpreted 
as the probability that the system has not emitted any 
photons since the previous emission event. 

In the event of a detector click in the or £)_ de- 
tectors, the state evolves discontinuously according to 

R+\4>) 



W+) 



\rIr+\^) ' 

R-\^) 



(15) 



(^|i?Li?-|V) ' 

respectively. The corresponding probability density for a 
click in either detector is given by 

u{t,±) = m)\RiR±m))- (16) 

Here, u{t^ zL)dt is the total probability that a click occurs 
in detector D± between times t and t + dt, and that no 
click occurred before time t. 

Following the first optical tt— pulse, the state of the 
two-atom system at t = is given by \ip{0)) = |(|00) -I- 
|02) -I- 1 20) -I- 1 22)). Note that we are interested only in 
the parts of the state that will ultimately be post-selected 
conditional on observing detector clicks on both rounds 
of the protocol. Thus we can restrict attention to the 
components of the state in the subspace spanned by the 
states |01), |10), |02), and |20). This leads to the fol- 
lowing closed set of coupled equations for the evolution 
generated by Eq. (fT4| . 



■A 

-zAco2 2~'^°^ 



C02 

C20 — 2~^^" ' 

coi = 0, 
Clo = 0, 

(17) 

where Cjk{t) — {jk\ip{t)). Eqs. pT)) may be readily solved 
to give 



C02{t) 

C2o(i) 



co2(0)e-('^+«"«/2)t 
C2o(0)e-''^'^*/2 ^ 



(18) 
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Since co2:oo(0) = C20;oo(0), we see that these two coeffi- 
cients are identical up to a varying phase factor, 



co2 = e ''^'c2o. 



(19) 



By applying Eq. (|15p . we find that the normalized zero- 
excitation component of the state after the first click at 
ti is given by 



|01)±e'^^^|10) 
V2 



(20) 



corresponding to a click in the D+ or D_ detector re- 
spectively. Note that the true state will include terms 
proportional to |12) and |21). However, these will be 
post-selected away on the second round of the protocol, 
and so can be neglected for the purposes of this analysis. 
Proceeding with the second round as described in Sec- 
tion |TT1 we find that the second click at t2 removes these 
unwanted terms, and leaves the system in 



mi,t2)) 



|01> + (-l)'V'^(*i-*2)|10) 

71 



(21) 



where m is the number of Z)_ clicks that have been ob- 
served. We note that regardless of the actual values of ti 
and t2, the final state for any detector click combination 
is always a maximally entangled state. In addition, as 
long as the values ti and t2 are known, it is possible in 
principle to undo this additional phase shift in the sys- 
tem with a local operation. This may be achieved via 
a rotation about the z-axis by an angle — A(ti — t2) on 
the first qubit. Therefore in this ideal case the fidelity 
of the protocol is unity despite mismatch of the cavity 
frequencies. 

This result can be understood as follows. In general we 
might expect a reduction in fidelity due to the fact that 
the frequency of the photons carries some 'which path' in- 
formation about which of the atoms the detected photons 
originated from. However, since we have assumed ideal- 
ized time-resolution detectors, complementary informa- 
tion about the frequency of the photon cannot be deter- 
mined, even in principle. Thus these idealized detectors 
themselves erase the 'which path' information, contained 
in the frequencies of the photons, which might otherwise 
have reduced the fidelity. 



B. Bad detector limit 

We now consider the case of very bad time resolution 
detectors, when tr ^ A~^, K~g. In this regime, the de- 
tectors give essentially no information about the arrival 
time of the photons, but simply indicate whether a pho- 
ton was observed or not on a given round of the protocol. 
We again assume that the collection and detection effi- 
ciencies are unity. 

By inspecting Eq. (PT]) it is clear that, since ti and 
t2 will now be unknown, an unknown phase will be ac- 
cumulated between the two terms in the superposition. 



Averaging over this phase will lead to a mixed state, with 
less than ideal fidelity. Assuming the click is observed in 
detector the state at the end of the first round of 
the protocol will be 



1 



P+ 



Pi+) Jo 
1 



dtip{ti,+) 



x|(|10)+e-''^'i|01))((10|+e'^*i(01|). (22) 

Here, p{ti, +)dt is the probability that a single photon is 
emitted into I?+ between times ti and ti + dt, and that 
no photon is observed subsequently in either detector. 
p(-l-) — Jg°° dtip{ti, +) is the total probability of observ- 
ing precisely one photon in D^, and no photons in Z)_. 
These quantities can be calculated within the QJ method 
to give 



P(ii,+) 



Kofie" 



-Kctftl 



(23) 



from which we obtain p{+) = 1/4. The state of the 
system at the end of the first round is found to be 



1 



P+ 



|01)(01| + |10)(10| 



KcS + «A 



|01)(10| 



Koff - «A 



|10)(01| . (24) 



To determine the effect of the second round of the proto- 
col, we use p+ as the input, and evolve under the gener- 
alized Schrodinger equation, p{t2) = — ■^[-ffcondP(^2) ^ 
p(i2)-ffcond]- The state corresponding to a click in 
D+ at time t2 is then i?+p+(t2)^+/tr[i?+/3+(f2)-R+], 
and the corresponding probability density is p{t2,+) = 

tr[i?+yO+(t2)^+]- From this we find that the final state, 
after both rounds of the protocol, is given by 



P++ = 



dt2R+p+{t2)Rl, 



i(|01)(01| + |10)(10| 



^eff 



:(|01)(10| + |10)(01|) . (25) 



A useful figure of merit for evaluating the effect of fre- 
quency mismatch is the fidelity of the final state with the 
closest maximally entangled state. Note that the proto- 
col is such that the final state only has components in 
the {|01), 1 10)} basis, so the closest maximally entangled 
state will be of the form \ip{<j>)) = (|01) -I- e''^|10))/\/2. 
Thus, for a given state p, the fidelity can be found by 
maximizing 



F(p,0) = (7/^(0)|p|V'(0)) 



(26) 



with respect to 
so we find that (, 
given by 



i.e. solving dF{p.(f)) / d(t> — 0. Doing 
= — arg(poi,io), and that the fidelity is 



np) 



iPoi.iol 



(27) 
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FIG. 2: Log-log plot of the final state fidelity as a function 
of Koff/A. We notice that even for a bad detector which is in 
principle able to distinguish the different frequency photons, 
a very high fidelity is possible provided that the bandwidth 
of the photons is large compared with their detuning. 

Substituting p++ from Eq. (|25|) gives 

We plot this result in Figure [H F(p++) increases toward 
unity for Keff ^ A. This result can be understood by 
considering the spectrum of the emitted photons. More 
specifically, Keff can be thought of as the bandwidth of 
each photon in the frequency domain. Thus for Keff 3> 
A, the photons become spectrally indistinguishable, and 
therefore carry no which-path information. In this regime 
we can obtain a large fidelity, even with very low time- 
resolution detectors. 



C. Intermediate detector case 

In this section we describe the effect of having a detec- 
tor with arbitrary time resolution, which may lie between 
the two limiting cases already studied in the preceding 
sections. This requires a more sophisticated model of the 
detection process to include a finite time resolution. 

The detector model we use is similar to the one de- 
scribed in Refs. 0,1111, ^^'^ summarized in Figure [S] 
While the complete model in Refs. incorporates 
a number of realistic imperfections, including finite time 
resolution, dead time, and dark counts, here we restrict 
our attention to the effect of finite time resolution alone. 

We now briefly review the model of [l^, HI] for the 
case of a single fluorescing quantum system, observed by 
a single detector, neglecting dead time and dark counts. 
In this case [l^, [11], the model consists of an idealized 
detector coupled to the transition between two internal 
states, ready and triggered (denoted here by R and T, re- 
spectively), of the detector. When a photon is registered 
by the ideal detector, the internal state changes from R 



Ideal Triggered, T 




Ready, R 



FIG. 3: Diagram showing the detector model adopted here to 
take account of the finite time resolution of the detector. 

to T. The T state of the detector decays, via a Poisson 
process, back to i? at a rate 7^. This decay process can 
be understood as the observed detector 'click', i.e. the 
time of the T ^ R transition corresponds to the pho- 
ton detection time reported to the observer. 7"^ can be 
understood as the response time of the detector. 

The quantities of interest are the unnormalized density 
matrices puit) and pxit). Pnit) corresponds to the state 
of the quantum system at time i, given that the detector 
is in the 'ready' state, while pt(0 corresponds to the 
state of the quantum system given that the detector is 
in the 'triggered' state. The states are unnormalized, 
such that tr[pfl(i) -t- prit)] is the total probability that 
no click has been observed up to time t (provided the 
state was normalized at time t = 0). The normalized 
state of the system at arbitrary times is then given by 

pit) = [PR{t) + PT{t)]/tT[pR{t) + pT{t)]. 

Between detector clicks, the state evolves according to 
the generalized mater equation 

PB. = J[R])PR^ , 

PT = iC~-/r)pT + J[R]pR, (29) 

where Cp = -■|(i?cond/0 - P-ffiond) + ^P^^ usual 
Lindblad superoperator, and J[R]p = RpR^ is the jump 
superoperator. R here is the appropriate lowering opera- 
tor for the fluorescing system of interest, with -ffcond the 
corresponding conditional, non-Hermitian Hamiltonian. 

The first equation corresponds to the conditional evo- 
lution for the case when no photon is emitted, and is a 
generalization of Eq. (fT4|) to mixed states. The second 
equation has three contributions. The term Cpx corre- 
sponds to the unconditional evolution of the system, i.e. 
it is what would be expected if there were no detector 
present. This can be understood since, once the detector 
is in the triggered state, no more information about the 
state of the system can be obtained from the detector un- 
til it returns to the 'ready' state via a 'click' event. The 
term proportional to "frPT ensures that pt is appropri- 
ately normalized, given that no click has been observed. 
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Finally, the J[R]pii term couples the two equations and 
corresponds to i? to T transitions of the detector when a 
photon is received. 

The probability distribution for observing a detector 
click at time t, given that no clicks were observed until 
time t, and provided the state was normalized at time 
t = 0, is given by 

=7rtr[pT(i)]. (30) 
When a click occurs, the state evolves discontinuously as 

PRit + dt) =PT{t), 

PT{t + dt)=0. (31) 

Eqs. ((29l - [3T|l allow one to calculate the probability distri- 
butions for click times, and also to determine the state of 
the system at arbitrary times, conditional on the record 
of detector click times. 

It is straightforward to generalize this model to the 
case of two detectors, observing two atoms via a beam 
splitter. We now have to track four quantities, prr, 
p+, P-, and ptt- Prr corresponds to the state of the 
two atom system when both detectors are in the ready 
state. p+ (p„) corresponds to the detector in the 
triggered (ready) state, and D- in the ready (triggered) 
state. PTT corresponds to both detectors being in the 
triggered state. 

The set of equations governing the evolution of the 
system conditioned on no detector click is 

PRR = {C-J[R+]-J[R-])pRR, 

P+ = {C- J[R-]--ir)p+ + J[R+\PRR, 
P- = {C- J[R+]--ir)p^ + J[R-]PRM, 
PTT = {C -2^r) P+-+J[R-]P+ + J[R+\P-. (32) 

I 



Poi,oi(ii, ^2) = 
^01,10(^1,^2) = 
Pio,io(iij ^2) — 



where ti and t2 are the times of the observed detec- 
tor clicks on the corresponding rounds of the protocol, 
and m = if both clicks are in the same detector, else 
m — 1. All other components of the density matrix van- 
ish. Note that, as before, this state is unnormalized, with 
the normalization such that the joint probability distri- 
bution for any specific two click combination is given by 
P{ti,t2) = 7rtr[p(ti, t2)]. As expected, this expression 
integrates to ^ which is the total probability of having 



where Cp = -j-{H,ondP - pfl"Ld) + R+pK + R-pRl 
is the Lindblad master equation, J[R]p = RpR^ is the 
jump matrix, and jr is the stochastic decay rate at which 
a triggered detector will signal a click and return into the 
ready state. 

The complimentary (instantaneous) evolution corre- 
sponding to a ZJ-i- click event is now 

PRRit + dt) = p+{t) , 
p+{t + dt) =0, 

p^{t + dt)^PTT{t), 

PTT(i)-0. (33) 
while a D_ click event is described by 
PRR{t + dt) = p_(t) , 

p+{t + dt)= pTTit), 

p-{t + dt) =0, 

PTT{t)=0. (34) 

We may solve the coupled equations of ([5^ analytically 
in the over damped regime, where -ffcond and R± are 
given by Eqs. ((T^] - [131) ■ The analysis can be consid- 
erably simplified by noting that we are only interested 
in the case where a single detector click is observed on 
each round of the protocol, since other outcomes will be 
discarded as failures anyway. This means that we need 
only to explicitly keep track of terms in the state which 
contain a single excitation. 

Evaluating the Eqs. (|52HM)) . for the case when exactly 
one detector click is observed on each round of the pro- 
tocol, leads to the final state of the system given by 



(35) 

I 

exactly two clicks in a specific sequence such as fol- 
lowed by 

For A ^ /toff, 7r we see that the coherence van- 
ishes corresponding to a large amount of mixing, which 
confirms the results from Section IIII Bl Similarly, for 
A <C KeS, Ir, we find that we recover the pure, maxi- 
mally entangled state from Section UlI Al 

More specifically, using Eq. ([77)1 we find that the fi- 
delity of the protocol for this intermediate detector case 



P20.20;o(0)Keff ,^-7r(ti-Ht2)^^~(Kcff-7r)ti _ i-^^^^-l,i^at{-l,)t2 _ 2) ^ 
4(Kcfr - 7r)^ 

4(iA -f Koff - IrY 

P02,02;o(0)Kgff ^-y.,(ti+t2)^^.-(K,ff~7r)tl _ 2)(e-('''"'f-'^')'=^ - 1) , 
4(Kofr ~ IrY 
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is given by 



0.5 



-(iA+K-7r)ti 



l)(e- 



-(iA + K-7r)t2 



1)1 



(k - 7,.)2 + A2 (e-(«-7.)ti _ i)(e-('--7.0t2 _ i) 



(36) 
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FIG. 4: (Colour online) Average final fidelity for the interme- 
diate detector regime, from a numerical evaluation of Eq. (I38[) . 



The average fidelity of the protocol is then 

- _ dh dhPihM) ■ F{hM) 



(37) 



Using Eqs. ([55)1 and (^7)) and the joint probability density 
for the clicks, P(ti,t2), this simplifies to 

F ^ h+ dti / dt287'|poi.io(ii,i2)|. (38) 



This expression is not readily integrable analytically, but 
it may be solved numerically. Doing so for a range of 
7r and KeS gives the results plotted in Figures H] and 
[5l Figure [4] indicates that the error in the entangling 
operation [that is, 1 — F, with F the average fidelity of 
Eq. p8|) ] scales approximately quadratically with A/jr, 
and also indicates that for fidelities exceeding 0.99 we 
require jr ^ 6 A. We also find that HeS and 7^ have 
similar influences on the final state fidelity, as can be 
seen from the plot shown in Figure [5l 



IV. HONG-OU-MANDEL DIP WITH DETUNED 
SOURCES 

In this section, we consider another few-photon inter- 
ference phenomenon, the Hong-Ou-Mandel (HOM) ef- 
fect, with detuned photons. The HOM effect is well un- 
derstood and was first observed in 1987 [6|. In the usual 



FIG. 5: (Colour online) Average final fidelity for the interme- 
diate detector regime, from a numerical evaluation of Eq. (|38[). 
as a function of both n^s and 7^. 



setup, two photons with identical spectra and spatio- 
temporal mode shapes simultaneously impinge on the 
two different input ports of a beam splitter. Provided 
these photons are identical and phase coherent, the state 
at the output ports of the beam splitter is then given by 
l^/j^m) = (|2,0) + |0, 2))/v^. In this idealized situation, 
subsequent detection of the photons by photodetectors 
placed at each output port thus reveals perfect bunch- 
ing or coalescence of the photons; that is, both detection 
events occur in the same detector, and coincidences cor- 
responding to one photon leaving each output port of the 
beam splitter are not detected. 

The HOM effect was originally used to characterize the 
coherence of individual photons from parametric down 
conversion sources @ , and has also been used to charac- 
terize photons from a single quantum dot source [40| . 
More recently, this effect has been used as the basis 
of many schemes for linear optical quantum computing 
[3, and has been proposed as a method for entan- 
gling spatially separated atomic ensembles The ef- 
fect has recently been observed using photons from inde- 
pendent sources, such as independently trapped neutral 
atoms [41'! and separate ions [3^. 

The ubiquity and utility of this phenomenon in 
quantum information processing therefore motivates the 
quantitative study of the HOM effect when the photons 
are not identical, in particular when the center frequency 
of the two photons is not identical. HOM interference of 
such detuned photons has already been studied both the- 
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oretically [i| and in an elegant experiment using photons 
from an atom-cavity system [sj. Legero et al. predicted, 
and subsequently observed, that perfect photon coales- 
cence is not seen. Rather, provided sufficient time reso- 
lution is available in the detector signals, the probability 
of detecting both photons in the same detector oscillates 
as a function of the detection time, with the frequency 
of this beat oscillation given by the relative detuning of 
the two input photons. 

Here, we study this effect using a similar method to 
that described in Section Hill i.e. we treat the effect as a 
continuous measurement problem, and keep track of the 
internal quantum state of the sources. In Section IIV Al 
we consider the limit of ideal detectors, namely ones for 
which precise information about the timing of the de- 
tection events is available. As noted above, a similar 
calculation has already been performed in [J using dif- 
ferent techniques. We obtain results in agreement with 
Ref. [1] , which demonstrates the validity of our approach, 
and furthermore, we hope sheds an alternative perspec- 
tive on the result. Then in Section flV Bl we discuss the 
opposite limit, and determine the coalescence probabil- 
ity in the case where no timing information is available 
from the detectors. We quantitatively consider the inter- 
mediate case in Section HV CI where the detectors have 
finite time resolution, and evaluate the visibility of the 
HOM beat fringes as a function of the both the photon 
detuning and detector resolution. 



A. Ideal detector case 

In this section we consider the HOM effect with de- 
tuned photons, of center frequencies ui and = t^i — A 
in the limit of detectors with time resolution much 
shorter than A~^. We use a simple model of the single 
photon sources, comprising single mode leaky cavities, 
each initially prepared in the single photon Fock state, 
as shown in Figure [6l 

As elsewhere in this paper, we treat this setup as a 
continuous measurement experiment, and apply the QJ 
formalism. We assume, for simplicity, that the efficien- 
cies for photodetection, collection, and photon emission 
into the desired mode are all unity. This assumption can 
be relaxed later. Provided the emission efficiencies are 
identical, the results will be essentially the same, except 
for an overall reduction in coincidence count rate. 

The quantities of interest are the time-resolved coin- 
cidence probability density functions, p{ti, t2, ±, ±). For 
instance, p{ti,t2,-\-,—)dtidt2 is the probability that the 
first click is observed in the detector _D_|_ in the infinites- 
imal interval [ti,ti -\- dti], and that the second click is 
observed in £)_ in the interval [^2, ^2 + dt2]- We restrict 
our attention here to the quantity p(ti,t2, +, +), but the 
other quantities are straightforward to calculate in a sim- 
ilar manner (in particular the symmetry of the setup im- 
plies p{ti,t2, +, +) = p{ti,t2, -, -)). 

Between photodetection events, i.e. conditional on no 



cavity 1 , 




cavity 2, oj. 



FIG. 6: (Colour online) Set-up for a HOM experiment with 
detuned photons. The single photon sources are each modeled 
as single mode cavities, with a leaky mirror of decay rate 
K at one end. The frequencies of the cavities (and hence 
the emitted photons) are uji and W2 = a;i — A respectively. 
The sources are initially prepared in the single photon Fock 
state |1, 1) = cjcj 1 0,0). Photons emitted from these sources 
impinge on the two different input ports of a beam splitter, 
and are subsequently detected by the photodetectors D+ and 
D-. The cavity sources may be replaced with a variety of 
different physical systems, with qualitatively similar results, 
as described in the text. 



detector clicks, the cavity systems evolve under the 
Schrodinger equation [Eq. ([M]) ] according to the non- 
Hermitian conditional Hamiltonian 



tiwib\bi 



ifiK i ihn 4. 
—b\b, - —blb2 . (39) 



Here, ujj is the frequency of cavity j (and hence the cen- 
tral frequency of the emitted photon) , bj are the corre- 
sponding lowering operators, and k is the leakage rate of 
each cavity. For simplicity, we assume that the cavity 
leakage rate, and hence the temporal mode shape of the 
emitted photons, is the same for each cavity. 

As in Section imi in the event of a detection event the 
system evolves discontinuously according to 

(40) 



(V|i?Li?-|v> 



(41) 



In this case, the jump operators are given directly by the 
beam splitter transformation. 



R+ = ^/ -(61+62). 



R- = \l^ibi-b2). 



(42) 



It is worth noting that while Eqs. ((39)) and (|42p . to- 
gether with the initial condition that each cavity is pre- 
pared in a Fock state, represent a rather idealized model 
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of a single photon source, they are also directly relevant 
for more realistic systems that could be used as sources. 
These include isolated two-level atomic systems (such as 
trapped atoms, ions, quantum dots, impurity centers in 
diamond), or corresponding systems coupled to a single 
(adiabatically eliminated) cavity mode in the bad-cavity 
regime. In such cases, the cavity mode is replaced with 
the corresponding atomic transition, and the lowering op- 
erators bj should be replaced with the corresponding low- 
ering operator for the atomic system, similarly to Section 
mil The results below therefore also apply to these alter- 
native systems. 

We assume the cavities are initially prepared in the 
Fock state |V'(0)) = |1,1) = b\bl\0,0), where |0,0) 
represents the vacuum mode for both cavities. Sub- 
sequently, conditioned on no detector clicks being ob- 
served between times and ti, the state for the cavities 



cond^OIV^O)) 



evolves according to |V'(^i)) ■ 
exp[— (iwi -I- iuj2 + K)ti]\l, 1). 

Assuming a click occurs at time ti in detector the 
state evolves according to Eg. 140]) . and we find that the 
normalized state after the detector click is 

rnti+dt)) = -L(b{ + bi) io,o). 



V2 

1 

71 



(|1,0) + |0,1)) . 



(43) 



The probability density for this detector click to occur 
between times ti and ti + dt, conditional on no clicks 
being observed up to time ii, is given by p{ti,+)dt — 
\\R+\'ipiti))\\'^dt = Ke-2''*idt. Note that Eq. ^ implies 
that the sources become entangled after the first detector 
click. Furthermore, the state of the system immediately 
after this first click is independent of both the click time 
ti and the detuning of the two photons, A. 

Conditional on no clicks being observed between times 
ti + dt and ^2, the state again evolves under the 
Schrodinger equation with iJcond, as 

\i;(t2)) -exp[-|i/eond(t2 - h)]\ij{h + dt)) , 

= ^ (e-(''^i + t)('2-ti)|l,0) 



-(ic^i + f )(t2-ti) 



V2 



_^g-(ic^2 + f)(t2-tl)|0^l)^ , 

(|1,0) +e'^(*^"*i)|0,l) 



(44) 



Note that the relative phase between the two terms in 
this expression oscillates in time at a rate A. It is this 
oscillation that leads to the 'quantum beats' in the pho- 
ton coalescence probability. 

A little algebra shows that the probability density for 
the second detector click to occur in detector be- 
tween times t2 and t2 -t- dt^ conditional on the first click 



A = 4n- 
A = 0.5tt 




FIG. 7: (Colour online) Generalized Hong-Ou-Mandel effect 
with detuned photons, (a) Probability density function for 
observing both clicks in the same detector, p{t2,+\tl,-\-) = 
p{t2, ~\tl, —), as a function of the delay between clicks, 
{t2 — t\). This probability density oscillates at a rate A, the 
detuning between the center frequencies of the two photons. 
K = 1, A = 47r (solid curve), A = O.Stt (broken curve), (b) 
Total coalescence probability, p(-l-|+) = p( — | — ) as a function 
of detuning, A/k. This corresponds to the coincidence signal 
that can be measured if no timing information is available 
from the detectors. 



being observed in at time ti , is given by 

p(i2,+|il,+)= ||i?+|V^(t2))|P, 

= ne 



"(*^-*^)-^{l + cos[A(t2-ti)]} 



(45) 



Combining this with the expression for p(ti,+) 
given above we find that the total joint probabil- 
ity distribution for two clicks in detector is 
given by p{ti,t2,+,+) = _p(i2, +1^1, +)p(ii, +) = 
^2g-K(t2-Hti) _^ cos[A(t2 - ti)]} /2. The somewhat 
more more instructive expression p(i2,+|il,+) is plot- 
ted as a function of t2 — ti in Figure [7] (a). 

The central result of this section is that the probabil- 
ity of finding both photons in the same detector oscillates 
with frequency A, the detuning of the two photons. This 
is in agreement with corresponding theoretical and ex- 
perimental results in Refs. [J| and [5[, although we have 
used a different method to arrive at the result |45|]. 
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B. Bad detector limit 

From the results of the previous section, it is straight- 
forward to calculate the total photon coalescence prob- 
ability - that is, the total probability that both photons 
will arrive in the same detector. This can also be thought 
of as the visibility of the Hong-Ou-Mandel effect in the 
limit of bad photodetectors, i.e. detectors whose time res- 
olution is much longer than the inverse of the detuning, 
A~^. We assume that in other respects, the detectors 
are ideal, in particular that their detection efficiency is 
unity, and no photons are missed due to dead time. In 
this limit, the only information available from the de- 
tectors is the total number of photons arriving at each 
detector, and no information about the time of arrival is 
available. 

The total photon coalescence probability, p{+\+) = 
p{—\—) for observing both clicks in detector is given 
by integrating Eq. (|45|) over all values of t — t2 — ti, 



P{+\- 



/ dTp{t2,+\tl, 

Jo 



dTKC '^^ • - {1 + cos At} 



A2 



(46) 



This expression is plotted in Figure [7fb). For small 
detuning, A ^ k, p(-|-|-f) approaches unity, indicating 
near perfect Hong-Ou-Mandel photon coalescence. For 
A ^ K, p(-f|-f) — 1/2, indicating that the photons 
apparently act as independent particles, scattering ran- 
domly into both output ports of the beam splitter, k can 
be thought of as the bandwidth of the photons, giving 
the spread in frequency of each photon about its center 
value. Thus these results agree with the usual intuition 
that HOM interference can not be observed between dis- 
tinguishable photons. However, it is important to note 
that this is a consequence of the poor time resolution of 
the detectors. As we saw in the Section flV Al with suf- 
ficiently fast time resolved detectors, a modified form of 
the HOM interference can be observed. 



Intermediate detector case 



In Sections HV^ and [TVBl we discussed the HOM ef- 
fect in the limiting cases of perfect time resolution pho- 
todetectors, and the limit of bad detectors (with no time 
resolution), respectively. In general, experiments may 
be in neither limit, and the detectors may have a finite 
time resolution which is comparable to the inverse of the 
detuning between the incident photons. In view of the 
applications of the HOM effect in quantum information 
processing, it is useful to have a quantitative understand- 
ing of the effect of finite time resolution detectors. In this 
section, we apply the model of finite time resolution de- 



tection introduced in llll CI to analyse the HOM effect in 
this regime. 

As in Scction llll Cl we model the detector by the set of 
equations 

PBR = iC-J[R+]-J[R-])pRB., 
p+ = {C-J[R4-^r)p+ + J[R+]pRR, 

p_ = (£- J[i?+]-7r)P-+ J[i?-]pflfl, 
PTT = {C-2-fr)pTT + J[R-]p+ + J[R+\p-. (47) 

where, as before, Cp = -\{Hcondp- pHl^^^) + R+pR\ + 

R-^pR'L is the Lindblad master equation, J[R]p = RpK^ 
is the jump matrix, and 7^ is the stochastic decay rate at 
which a triggered detector will signal a click and return 
into the ready state. Now, however, the system under 
observation is that shown in FigureEl with -ffcond is given 
by Eq. ([39)1 , and the operators R± are given by Eq. (I42p . 
Thus the matrices prr^ p±, ptt describe the internal 
states of the two single mode cavities representing the 
sources. As in Section flV Al we take the initial state of 
the system to be the Fock state |1, 1), such that prr{0) = 
|1, 1)(1, 1|, and p±{0) = ptt{0) = 0, indicating that both 
detectors are in the 'Ready' (R) state at time t = 0. 

We are interested in quantities such as p^(ii, +, +), 
which is the probability density function for pairs of ob- 
served detector clicks. This can be found by first solving 
Eqs. ([Tf|) to find Pj{ti). At a given time ti, the total 
probability for detector to be in the 'Triggered' (T) 
state is Ts:[p+{ti) + pTxiti)]- Thus the probability den- 
sity for the first detector click is given by 

P^{ti,+) ^ JrTr[p+{ti) + pTTih)] . (48) 

After the first detector click, the state of the sys- 
tem is updated and renormalized according to 

p+{h) 



PRRih + dt) = — 

Tr[p+(ti) + pTT\ti)\ 

p+{ti+dt) = 0, 

P-[ti + dt) = — — TTTi ' 

TT[p+[ti) + PTT[tl)\ 

prriti + dt) = . (49) 

This state is then used as an initial condition to Eqs. l47l 
which can be solved to give pj{t2\ti, +), i.e. the state of 
the system at time t2, conditional on a £>+ click at time 
ii, and no clicks between times ti and ^2- The conditional 
probability density for the second click, in detector D+, 
is then given by 

P-y{t2, +\tl,+) = 7^Tr[p+(t2|il, +) + PTT{t2\tl, +)] . 

(50) 

The joint probability distribution for _D+ clicks at time 
ti and t2 can then be formed from Eqs. and (|50p as 
P-y{ti,t2, +, +) = P-y{t2, + 1*1 > +)P7 (^1 , +) ■ Notc th at, un- 
like in the ideal detector case described in Section IIV Al 
P^{t2, +1*1, +) is no longer just a function of ^2^*15 but 
in fact depends on the value of ti. 
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In principle, quantities such as p^(ti,t2, +, +) can be 
obtained analytically by first finding the general solution 
of Eqs. ([47]) and then applying the above steps. How- 
ever, this turns out to be difficult in general. Unlike 
the case of the distant-atom entangling scheme presented 
in Section fill B I it is not possible to concentrate on the 
single photon terms in the expressions for pj. Instead, 
all terms must be tracked, including in particular those 
which correspond to both photons entering the same de- 
tector within a short interval. This means that all 64 
coupled equations in Eqs. (|T7|) must be solved. Owing 
to this complexity, we resort to a numerical solution of 
Eqs. (|47)) . using the Matlab 'ode45' function. Follow- 
ing the steps outlined above, we obtain i2, +, +)• 
From this quantity, we can determine (r, that is, 
the distribution of intervals t = t2 — ti between detector 
clicks, given that the first click was observed in detector 
£)_!_. This is the analogous quantity to ^(^2, +1^1, +), in 
the case of ideal detectors, which was plotted in Figure 
□ (a). p^{t,+\+) is given by 



dtip^{ti,ti + r, +, +) 



(51) 



We can also calculate Pj{t, — \+), which is the proba- 
bility distribution for the second click to occur in Z?_, 
given that the first was in _D_|_, in a similar way. In prac- 
tice, we set the upper limit of the integrals to a value 
of 3k~^, which leads to a slight relative underestimate 
of approximately exp[— 2k x 3k^^] « 0.25% in the value 
of these integrals, ^^(t, -|-|-|-) is plotted in Figure Elja), 
along with the ideal result, p{t2, +\tl, +) from Eq. 
for the same values of A and k. For smaller values of 
the detector response rate, 7r, it can be seen that the in- 
terference fringes become washed out. This reduction in 
visibility is a result of the finite time response of the de- 
tectors. More specifically, as jr approaches A it becomes 
harder to resolve the the individual interference fringes. 

A reasonable figure of merit for this generalized HOM 
effect is the fringe visibility. 



v{t) 



(52) 



p^(T,+|+)+p^(r,-K 

After an initial transient behaviour on short timescales 
T ~ Ir^, v(t) is, to a good approximation, a peri- 
odic function, with maxima at integer multiples of tt/A. 
A good figure of merit is therefore the value of these 
maxima (taken outside the initial transient region). In 
Figure [5Kb) we plot the fringe visibility evaluated at 
T = 20 X tt/A, as a function of detector response rate, 
7r, for a particular choice of A and k. As expected, 
the fringe visibility vanishes for slow detectors, and ap- 
proaches unity for fast detectors with 7^ ^ A. 



V. OUTLOOK AND CONCLUSIONS 

In this paper we have looked at the effect of frequency 
mismatch in two schemes which make use of few-photon 
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FIG. 8: (Colour online) Generalized Hong-Ou-Mandel effect 
with detuned photons and finite time response detectors, (a) 
Probability density function for observing both clicks in the 
same detector, p-y{T, -l-|-l-), as a function of the delay between 
clicks, r = t2 — ii. k = 0.1, A = 2it, 7^ = 2 (solid curve), 
7r = 10 (broken curve). The dotted curve shows the re- 
sult corresponding to ideal detectors, i.e. p{t2, +\tl, +) from 
Eq. (|45p . for the same values of k and A . (b) Fringe visibility 
at T = 20 X tt/A as a function of detector response rate, 7r, 
for the same k and A as in (a). 



interference. Our main conclusions can be summarized 
as follows. 

We first looked at the effect of frequency mismatch in 
the distributed QIP scheme of Barrett and Kok |lj] . One 
key finding was that, with idealized detectors (i.e. such 
that the detector response rate 7^ is much faster than 
the photon detuning A) perfect entanglement of distant 
atoms is still possible. The entangled states pick up a 
known phase that depends on the times of the photode- 
tection events, and which can be corrected via a local 
unitary operation. This is true as long as the cavity mis- 
match is known from, e.g. previous measurements on 
the cavities. Any unknown variations in this parameter 
cannot be rectified. In the opposite limit of poor time 
resolution detectors (7,. <C A) the fidelity of the result- 
ing mixed entangled state depends on the overlap of the 
photons' spectra according to Eq. 

With the aid of a model of the internal dynamics of 
the photodetector, we were also able to analyse the in- 
termediate case, where the detector response rate, 7^ is 
of similar order to A. Here, we found that as 7^ exceeds 
A, the fidelity of the entangled states approaches unity. 
The numerical results suggest that the error in the en- 
tangling operation (that is, 1 — -F, with F the average 
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fidelity of Eq. ([38)) ) scales approximately quadratically 
with A/7r, such that an error of 1 — ^ 10~^ is possible 
with A/7r - 10"2_ 

We also analysed the Hong-Ou-Mandel effect with de- 
tuned photons, and arrived at similar conclusions. In 
the case of idealized detectors with time resolution much 
larger than the photon detuning, a modified Hong-Ou- 
Mandel effect can be observed, with the coalescence prob- 
ability oscillating as a function of time, at angular fre- 
quency given by A, the detuning of the two sources. This 
result is in agreement with earlier theoretical and exper- 
imental results by Legero et al. 0, [1] , although our ap- 
proach was slightly different to this earlier work. We also 
extended the analysis to the case of non-ideal detectors, 
and found that the fringe visibility of this modified Hong- 
Ou-Mandel effect decreases dramatically as 7,- is reduced 
below A, but approaches unity for 7^ 3> A. 

The usual understanding of few-photon interference ef- 
fects is that they can only be observed with identical 
photons. This is because experiments such as those con- 
sidered here make use of a beam splitter to coherently 
erase 'which path' information. Thus any 'excess' infor- 
mation, such as frequency or polarization, carried by the 
photon, that can be used, even if only in principle, to 
infer the previous path of the photon, will typically de- 
grade the interference effect. Even if the detector has 
no output corresponding to the energy of the incident 
photon, we could imagine that this information might be 
encoded in the internal state of the detector, e.g. in the 
energy of the exciton-hole pair created when the photon 
was initially absorbed, and so it could be determined, in 
principle, with a careful measurement of the microscopic 
state of the detector. 

This is in apparent contradiction with the results of 
this paper, which predict that a modified version of such 
interference effects can be observed with high fidelity 
even when the photons are spectrally distinct. How can 
the results here be reconciled with the conventional un- 
derstanding described in the previous paragraph? One 
way of understanding this is that detectors with suf- 
ficiently good time resolution are themselves capable 
of erasing some of this 'excess' frequency information. 
Loosely speaking, for a detector with time resolution 7"^, 
a detected photon is localized in a time window of corre- 
sponding duration 7"^. Thus the frequency uncertainty 
of the detected photon must be equal to or larger than 
7r. This means that there is no way, even in principle, of 
determining the frequency of the detected photon. Thus 
the 'excess' frequency information has been successfully 



erased. 

Our results should also shed light on the question of 
where the interference 'goes' in few photon interference 
experiments with detuned photons. With idealized de- 
tectors with high temporal resolution, we have seen that 
a modified version of the interference or entanglement is 
present, and in this sense there is no loss of interference 
visibility or entanglement. With imperfect detectors, the 
corresponding effects are diminished. Therefore another 
way of understanding these results is to say that imper- 
fect detectors introduce noise into the interference pat- 
tern, because the timing of the detector clicks produces 
additional randomness above what would be expected 
from spontaneous emission alone. 

We conclude by suggesting some possibilities for future 
work on the subject of few photon interference with de- 
tuned photons. One immediate application of our work 
is that it can be used to place design constraints on near 
future experiments - given detectors with a particular 
temporal resolution (as quantified by 7^ in this work), we 
know that the sources must be tuned such that A < 7,. 
in order to see the corresponding interference or entan- 
glement effect with reasonable visibility/fidelity. Con- 
versely, if tuning the sources is difficult, we know the 
corresponding detector resolution required. 

Similarly, for scalable quantum computation, the re- 
sults here can be used to place constraints that the phys- 
ical parameters A and jr must satisfy in order to imple- 
ment a fault-tolerant computation scheme, either in a hy- 
brid matter-light setup |lj, or in a linear optical scheme 
0, H, [13 ■ this case, some further work is required 
in order to relate the physical errors due to frequency 
mismatch and finite detector response rate to the more 
abstract error models employed in the analysis of fault 
tolerant computation schemes. 

Finally, although we have concentrated on two partic- 
ular examples of few photon interference in this paper, 
the results presented here, the methods used to arrive at 
them, and the resulting physical insights, should carry 
over to many other few photon interference setups. We 
therefore hope these results have useful implications for 
quantum information processing schemes which utilize 
few photon interference. 
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